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We use a bosonization approach to calculate the single- 
particle Green's function G(r, r) of non-relativistic fermions 
coupled to transverse gauge-fields in arbitrary dimension d. 
We find that in d > 3 transverse gauge-fields do not destroy 
the Fermi liquid, although for d < 6 the quasi-particle damp- 
ing is anomalously large. For d — * 3 the quasi-particle residue 

t( — ) 2 ], where k is the Thomas- 



vanishes asZa expf— rr-r-, s — , 

r 1 2-K{d— 3) v mc ' 

Fermi wave-vector, m is the mass of the electrons, and c is 
the velocity of the gauge-particle. In d = 3 the system is a 
Luttinger liquid, with anomalous dimension 7_l = {^(rnc') 2 - 
For d < 3 we find that G(r, 0) decays exponentially at large 
distances. 
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More than 20 years ago Holstein, Norton and Pin- 
cus realized that the self-energy of non-relativistic 
fermions that are coupled to the transverse radiation field 
is plagued by logarithmic singularities. They concluded 
that in d = 3 dimensions such a system cannot be a 
Fermi liquid. However, due to the smallness of the finc- 
structure constant — ~ -h^f the deviations from conven- 
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tional Fermi-liquid behavior are of little practical conse- 
quences. The discovery of unusual normal-state proper- 
ties in the high-temperature superconductors Q and re- 
cent investigations of half-filled quantum Hall systems ||] 
have revived the interest in this problem. In these sys- 
tems the relevant dimensionless coupling constant may 
not be small, and the effective dimensionality is not nec- 
essarily d = 3. Thus we are lead to consider the problem 
of studying electrons that are minimally coupled to a 
gauge-field — [4>, A]. In an Euclidean functional inte- 
gral approach, the model is defined in terms of an action 
S = Si m + S 2 {V, A^} + S 3 {A^}, with 



Si {^} = (3 ^2 % [-iu> n + £k] i>k , 
k 



-q^q + h q 1 A- q - A q 



(1) 

(2) 
• (3) 

is the Euclidean action of non-relativistic spinless elec- 
trons, which are represented by a Grassmann-ficld with 
Fourier components tpk- The excitation energy relative 

i 2 

to the chemical potential /i is £k = 2m ~ Mi where m is the 
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mass of the electrons. The coupling between matter and 
gauge-field is described by S2{tp, A^}, where g is the rele- 
vant coupling constant. The composite Grassmann- fields 

p q = Ylk V'feV'fc+g an d i q represent density and gauge- 
invariant current density. Finally, 5*3 is the action 
for the gauge-field in Coulomb gauge, q- A q = 0. For the 
Maxwell-field the functions / q and h q are 



/q Vq 2 



7 = P 4TT 

9 I/q2 + (^)2 



(4) 



where V is the volume of the system, (3 is the inverse 
temperature, and c is the velocity of the gauge-particle. 
Throughout this work we shall use the convention that 
k = [k, iuj n ] and q = [q, iui m ] , where the fermionic fre- 
quencies are Co n = 2n(n+ \)/ P, and the bosonic ones are 
Lj m — 27rm//3. 

The integration over the transverse gauge-field gives 
rise to a retarded current-current interaction between the 
matter degrees of freedom |^,|[ . In d < 3 this interaction 
cannot be treated within the usual perturbative many- 
body theory, so that it is necessary to attack the problem 
by means of a non-perturbative approach, such as renor- 
malization group and scaling methods [^|- generalized 
jr expansions || , higher dimensional bosonization meth- 
ods p0| , or other infinite resummation schemes pl|Jl2] |. 

The bosonization technique in d > I is based on ideas 
due to Luther jl3l . It was rediscovered by Haldane H| , 
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and has been further developed in Refs. |]l5| , |l6 
recently, Schonhammer and the present author [1 
developed a powerful formulation of higher-dimensional 
bosonization which is based on functional integration. 
In the literature on one-dimensional bosonization these 
functional techniques have been known for some time 
p8[ , but have not received much attention. In Ref. |ll]] 
we used our approach to calculate corrections to the 
non-interacting boson approximation, which exist even 
in d = 1 if the energy dispersion is not linearized. Up 
until now these corrections could not be calculated by 
means of the conventional operator approach In 
the present work we shall show that the problem of elec- 
trons coupled to transverse gauge-fields can be treated 
within our method in a very simple and elegant way. 

We start by partitioning the Fermi surface into patches 
of area A d_1 and subdividing the degrees of freedom 
close to the Fermi surface into squat boxes K a of hight 
A associated with the patches |IJ,|T5|,|T^,^9| . The boxes 
and patches are labelled by an index a in some con- 
venient ordering. In Ref. we have decoupled the 
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two-body density-density interaction via a generalized 
Hubbard- Stratonowich field <f> q which couples to the lo- 
cal density-field p q = J^k ® a (k)V'!V'fc+gi where 0°(k) is 
unity if the wave- vector k lies inside the box K a , and van- 
ishes otherwise. Note that by construction p q = J2 a Pa,- 
Imagine now that in our gauge action defined in Eqs.|i|- 
U we integrate first over the 0-field, decompose then 
the total density into its contributions from the vari- 
ous patches and finally decouple again the matter action 
by means of the auxiliary field cf> q introduced in Refs. 
p7| , ^9[ . Clearly the longitudinal part of the new decou- 
pled action can be obtained from Eqs.|| and by replac- 
ing Y.qPq^-q ~> Eg Ea Pg4>-q and Eg/g ^-q^q ~> 

Eg HaALq \ aa '> where f_g is a matrix in the 
patch labels, with matrix elements [f q ] aa = / q [p0| . 
Hence, the Hubbard-Stratonowich field 4> q that was for- 
mally introduced in Ref. Jl7|] to decouple the two-body 
density interaction can be identified physically with the 
patch-version of the longitudinal gauge-field. Evidently 
the transverse part of the gauge-field can be treated in 
exactly the same way, so that the correlation functions of 
the "0-field can be obtained by functional integration over 
the patch gauge-field action S = S {ijj} + §2 {tp,A"} + 
§3 {Afy. The form of §2 and S3 is similar to Eqs.|| and 
||, the only difference being that the gauge-fields and the 
densities carry an extra patch-label a. Thus, the second 
terms in Eqs.|| and || are replaced by E 9 E Q jg ' A" 9 

and E 9 Eq Q '[^ ? 1 } aa ' AZ q ■ A™', where h q is a matrix 
in the patch labels, with matrix elements [h q ] aa = h q 
|§. The patch current density is j° = f q araa + jf a < a 

With % ara > a = £ fe e«(k)£±^)^fc+g and jdto,a = 
2mc 2 P i-iq' ^-q-q'Vq'- 

To calculate the fermionic Green's function G(k), we 
start from the representation of G(k) as a functional in- 
tegral over the fields ip and A^, and integrate first over 
the ip-tie\d. In this way we obtain the Green's function 
Q(k, fc'; {A^}) of the matter field for a frozen config- 
uration of the gauge-field. Note that this function is 
not diagonal in momentum-energy space, because the 
gauge fields act like a dynamical random-potential, which 
breaks the translational invariance in space and time. 
The full fermionic Green's function can then be obtained 
by averaging G(k, k] {A°}) over all configurations of the 
gauge fields, 



G(k) = (g(k,k;{A^})) j 



(5) 



where < . . . >a denotes functional average with the ef- 
fective gauge-field action S g f {Afy that contains in addi- 
tion to 5 3 {A™} a contribution Skin{A™} describing the 
influence of the matter degrees of freedom on the dy- 
namics of the gauge- field. Note that G(k, k;{A^}) = 

[G]kk, where the infinite matrix G is defined by writing 
SiM+W, A%} = -PEkk' ^{[G-^kk^k'- Of course, 
in general we cannot invert the matrix G -1 . However, if 



the Thomas- Fermi wave- vector k = {An g 2 is small 
compared with kp, then the matrix G" 1 is dominated 
by a narrow band of width k <C kp in the wave- vector 
label. Here v is the density of states at the Fermi energy, 
and kp is the Fermi wave-vector. If we then partition 
G -1 into block-matrices corresponding to the patches 
and choose the cutoffs A and A large compared with n, 
it is obvious that G _1 is approximately block-diagonal, 
[G~V' = E a a (k)e a (k')[(G«)-V', with diagonal 
blocks (G Q ) _1 associated with the boxes K a . Note that 
scattering processes which transfer momentum between 
different K a are neglected in this approximation, so that 
possible non-trivial effects due to the global topology of 
the Fermi surface are not taken into account |2lf| . 

The crucial point is now that within a given block we 



may linearize the energy dispersion, £k° 



v a -q, where 



k" is a vector on the Fermi surface that points to the suit- 
ably defined center of K a . Although each diagonal block 
(G")" 1 is still an infinite matrix in frequency space, it can 
then be inverted exactly by means of a method originally 
due to Schwinger p2fl . Thus, the matrix elements [G a ]tt 
can be explicitly calculated as functional of the A" . Fi- 
nally we average this expression over all configurations 
of the gauge-field, using the Gaussian approximation for 
the effective gauge-field action. As shown in Ref. [fl9|| , 
in a model with linearized energy dispersion the Gaus- 
sian approximation becomes very accurate at high densi- 
ties due to a large-scale cancellation between self-energy 
and vertex corrections (generalized closed loop theorem). 
The Coulomb gauge condition q • A q is easily imposed 
in our functional integral approach via the Fadeev-Popov 
method p3| . It should be kept in mind, however, that 
the single-particle Green's function is gauge-dependent, 
although physical observables derived from it are not. 
The final result for the interacting Green's function is 



G(/fc)=Ve Q (k) [ dr [ dre 
x G£(r,T)e Q » (r ' r)+Q ° (r ' r) 



— i[(k— k Q )-r — uj n r\ 



(6) 



Qf (r , r) = iq - Sf (r , r) , JJf = lim Q S|f (r , r) , (7) 
Ql(r,T)=Rf-Sl(r,T) , R a ± = lim Sl(r,r) , (8) 

r,r— »0 



[3V 



E 



cRPA.a 



cos(q • r - uj m r) 



(iuj m - v« • q) 5 



q<*(r —\ 1 V H * PA,a C ° S(q ' r Z ^"' T) 

PV 2 ? (*o; m -v«-q) 2 



G^(r,r)=4 d - 1) (r ± ) 



2w I r\\ + i\v a \T 



(9) 
(10) 

(11) 



where r« = r • v Q , with v" = v"/|v Q |. For length 
scales |rj_| ^> A -1 the function (5^ d ~^ (rjj can be treated 
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as a d — 1-dimensional Dirac-(5 function of the compo- 
nents of r orthogonal to v". Here f q HPA ' a is the effective 
screened interaction due to the longitudinal component 
of the gauge-field, as defined in Rcfs. |]l7| , ^9| . The term 
Q^(t,t) has been discussed extensively in our previous 
work fl7jj . Here we focus on the contribution Q" (r, r) 
due to the transverse part of the gauge field. In Ref. E4] 
we shall give an expression for h PPA ' a for an arbitrary 
shape of the Fermi surface and general patch-dependent 
bare interactions [Zi g ] QQ ■ For simplicity, in this discus- 
sion we shall restrict ourselves to a spherical Fermi sur- 
face (so that |v Q \ — vp for all a) and constant interaction 
matrices [^ g ] aa = hq, where h q is given in Eq.||. In this 
case 



R.PA,a _ 



1 



q) 2 



k) + ( V c ) ( VpK ) 2 + 9d ( 



(12) 



where q = q/|q|, and the function g d {z) is given by 



9d{z) 



-ld,z 



* d6 



7T COS 6* 



(13) 



Here j d 



Si,, 



g^j- n , where Cl d = ^ 2 /r(§) is the surface 
of the (i-dimensional unit sphere. Note that 7^ 



d — > 1. Because Gq(t,t) is proportional to 6^ ^(rj_), 
we may set r = rii v Q in the rest of the integrand of EqJq. 
Furthermore, close to and below three dimensions the 
value of the integral in Eq.|l^ is dominated by the regime 
C 1, so that we may approximate in Eq.|Tg| 



I for 



9d 



( i^>n 



\vfH 



Id- 



VF\q\ 



(14) 



Note that to this order it is consistent to neglect the 
term of order uj^ in the denominator of Eq, 12. After 
some trivial rescalings we obtain in the limit V, ft — ► 00 



1 



d-l 



dpp 



d-2 



COS 



du- 



d0(sm6) d \cos9\ — 
cos8\ipKd(r\\ —usovft) 



[p 2 + \u\] [lU - Sg] 



(15) 



where Cd = f^r (~fd) 2 is a numerical constant, Kd = 

1 

7? — k is the relevant wave- vector scale for the transverse 
gauge-field, and sg = sign{cos9). 

From now on we shall restrict ourselves to r = 0. The 
it-integration in Eq.O is then easily carried out, 



■2 



-c d 



d-l 



Tip" — p 



2 / dpp' 
mcJ J 

-l + (p 4 -l)lnp 2 



d-2 



[1+P 4 ] 2 



A d {pn d \ n \) , (16) 



where the dimensionless function Ad(x) is given by 

1 — cos ^| cos6>| 5 i 



A d (x) = / d6>(sin(9) d |cos6>| 
Jo 



(17) 



Obviously the large-r|| asymptotics of (r|| v a , 0) is de- 
termined by the behavior of A d (x) as \x\ — > 00. In d > 3 

For d < 6 the 

d-3 



the limit Af 



linix—joo A d (x) exists. 



leading correction is A d (x) ~ A c a 



with 



A?, 



r(^)W) 

r(M _i) 



d-3 



^=3^ 3 



,Trd. 
)sm(—), (18) 



It follows that for K d \r\\ 



oc 



Qi(niv Q ,o) 



-a d ' - 



, -) 

Kmc/ 



(19) 



with a,d — A^BdCd, where 



f°° d 9 irp 2 -p 4 -l + (p 4 -l)\np 2 , . 

s " = 2 i dw n+& ' (20) 

The integral is convergent for 1 < <i < 5 and can be 
evaluated analytically [Q . For d = 3 we obtain B3 = | • 
Since i?" is finite, transverse gauge-fields do not de- 
stroy the Fermi- liquid in d > 3: their contribution to the 



quasi-particle residue is Z? 



However, in d < 6 



the damping is anomalously large, since the leading cor- 
rection decays slower than | a; | 1 |17|] . In the limit <i — > 3, 



C 3 = -f- and A°, 



so that 



Z"=e R ° «e"M^) ( ™ )2 , 0<d-3«l . (21) 

In the marginal case d = 3 it is easy to show that 
A^{x) ~ I ln(|a;|) for |x| — + 00, so that for K3 |r|| | — > 00 

g«(r||V a ,0)^-|B3C 3 (— ) 2 ln( K3 |r|||) . (22) 
" 3 Kmc/ v 11 ' 



Exponentiating this expression yields 



eOKr^-.o) a (« 8 |r|,|)- 



•7x 



«3H|I 



■DC 



(23) 



where the anomalous dimension is 7 1 = J- f— ) 2 . In the 
case of electromagnetism at realistic densities — « 10 -2 , 
so that 7j_ ~ 5 x 10 -6 . Although this is extremely small, 
at distances large compared with k^ 1 the true asymp- 
totic behavior of the Green's function is characterized 
by a non-universal power law. This is precisely the Lut- 
tinger liquid behavior familiar from one-dimensional sys- 
tems with regular interactions. 

In d < 3 the function A d (x) diverges for \x\ — * 00 
stronger than logarithmic, A d {x) ~ ^l^l^ - with A' d — 
3 ^ 3 r(|)sin(^). It follows that for n d \r\\\ -> 00 



Qx(niv Q ,o) - -a d (— ) 



d-l 



Kmc J 



3 



(24) 
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with dd — A!' d B d+ 3Cd- Thus, for d < 3 the Green's 
function decays at equal times and large distances faster 
than any power. Obviously the quasi-particle picture is 
not valid but the non-universal power laws characteristic 
for conventional Luttingcr liquids do not appear. Such a 
system can be called a generalized Luttinger liquid. 

We would like to emphasize that our expression for 
<3"(r,T) in Eq.[l5| and the subsequent equations derived 
from it are independent of the box cutoffs A and A. 
This is due to the fact that k is the relevant ultravio- 
let cutoff for q-integrations. Because we have assumed 
k <C A, A <C kp, the integrals do not depend on the pre- 
cise choice of A and A. In this regime the bosonization ap- 
proach is most powerful and leads to cutoff-independent 
results that involve only physical quantities. On the 
other hand, for A,A < k the q-integrations extend over 
the entire volume of the boxes K a , so that the cutoff A 
appears explicitly in the final expression for the Green's 
function p0| . The marginality of d = 3 is also evi- 
dent in renormalization group calculations , but our 
bosonization approach goes beyond these works, because 
the calculation of the full Green's function in arbitrary d 
has been reduced to a well defined purely mathematical 
problem: After substituting EqJlJ into Eq|| the Fourier 
transformation can in principle be carried out numeri- 
cally, or perhaps even by some sophisticated analytical 
method. Finally, it should be mentioned that in Refs. 
P,[l2[ the validity of the bosonization approach in the 
gauge-field problem has been questioned. Note, however, 
that the generalized closed loop theorem jl9| guarantees 
that the cancellations between self-energy and vertex cor- 
rections, which are responsible for the accuracy of the 
Gaussian approximation at high densities, happen also 
in the case of gauge- fields Whether non-linearities 
in the energy dispersion modify the asymptotics of the 
Green's function is currently under investigation. We 
have preliminary evidence that this is not the case p4| ]. 

In summary, we calculated the Green's function of elec- 
trons coupled to the Maxwell-field at high densities. In 
the absence of spontaneous symmetry breaking the sys- 
tem is a Fermi liquid in d > 3, a Luttinger liquid in 
d = 3, and a generalized Luttinger liquid in d < 3. Our 
non-perturbative method is very general and can be used 
to study the coupling between electrons and any other 
bosonic field, such as phonons or magnons 

I am grateful to K. Schonhammer for many discussions, 
fruitful collaborations, and comments on the manuscript. 
I would also like to thank G. E. Castilla for detailed com- 
ments on the manuscript and W. Metzner for discussions. 



[2] 
[3] 

[4] 

[5] 
[6] 
[7] 
[8] 

[9] 

[10] 

[11] 

[12] 
[13] 
[14] 

[15] 
[16] 

HZI 

PI 
[19] 
[20] 



C. J. Pethick, G. Baym and H. Monien, Nucl. Phys. A 
498, 313c (1989). 

L. B. Ioffe and A. I. Larkin, Phys. Rev. B 39, 8988 
(1989); P. A. Lee and N. Nagaosa, Phys. Rev. B 46, 
5621 (1992); B. Blok and H. Monien, Phys. Rev. B 47, 
7312 (1993). 

B. I. Halperin, P. A. Lee and N. Read, Phys. Rev. B 47, 
7312 (1993). 

J. Gan and E. Wong, Phys. Rev. Lett. 71, 4226 (1993). 
J. Plochinski, Nucl. Phys. B 422, 617 (1994). 

C. Nayak and F. Wilczek, Nucl. Phys. B 417, 359 (1994). 
S. Chakravarty, R. E. Norton and O. F. Syljasen, Phys. 
Rev. Lett. 74, 1423 (1995). 

L. B. Ioffe, D. Lidsky, and B. L. Altshuler, Phys. Rev. 
Lett. 73, 472 (1994); B. L. Altshuler, L. B. Ioffe and A. 
J. Millis, Phys. Rev. B 50, 14048 (1994). 
H.-J. Kwon, A. Houghton, and J. B. Marston, Phys. Rev. 
Lett. 73, 284 (1994); preprint (Jan. 1995). 

D. V. Khveshchenko and P. C. E. Stamp, Phys. Rev. 
Lett. 71, 2118 (1993). 

C. Castellani and C. Di Castro, preprint (1994). 

A. Luther, Phys. Rev. B 19, 320 (1979). 

F. D. M. Haldane, Helv. Phys. Acta. 65, 152 (1992); 

in Perspectives in Many-Particle Physics, Proceedings 

of the International School of Physics "Enrico Fermi", 

Course 121, (North Holland, Amsterdam, 1994). 

A. Houghton et al, Phys. Rev. B 50, 1351 (1994); ibid. 

48, 7790 (1993); J. Phys. 6, 4909 (1994). 

A. H. Castro Neto and E. Fradkin, Phys. Rev. Lett. 72, 

1393 (1994); Phys. Rev. B 49, 10877 (1994); Phys. Rev. 

B 51 (in press). 

P. Kopietz and K. Schonhammer, preprint cond 
mat/9502088[ (July 1994). 



H. C. Fogedby, J. Phys. C 9, 3757 (1976); D. K. K. Lee 
and Y. Chen, J. Phys. A 21, 4155 (1988). 
P. Kopietz, J. Herm isson and K. Schonhammer, preprint 
cond-mat/9502089|, (February 1995). 



[21] 
[22] 
[23] 

[24] 



Of course, a matrix with constant elements is not invert- 
ible. However, the appearance of / and is only an 
intermediate step. Formally, we assume that we have reg- 
ularized / and hq in some convenient way and remove 
the regularization at the end of the calculation. 
S. Khlebnikov, Phys. Rev. B 50, 6954 (1994). 
J. Schwinger, Phys. Rev. 128, 2425 (1962). 
See, for example, G. Sterman, An Introduction to Quan- 
tum Field Theory, (Cambridge, 1993), p. 190. 
P. Kopietz, Habilitationsschrift, 1995, in preparation. 



[1] T. Holstein, R. E. Norton, and P. Pincus, Phys. Rev. B 
8, 2649 (1973). 



4 



